Efficient optical pumping of Zeeman spin levels in Nd^"^ : YVO4 

Mikael Afzelius*, Matthias U. Staudt, Hugues de Riedmatten, Nicolas Gisin 

Group of Applied Physics, University of Geneva, CH-1211 Geneva 4, Switzerland 

Olivier Guillot-Noel, Philippe Goldner, Robert Marino, Pierre Porcher 

Ecole Nationale Superieure de Chimie de Paris (ENSCP), Laboratoire de Chimie de la Matiere Condensee de Paris, CNRS-UMR 7574, ENSCP, 

11 rue Pierre et Marie Curie 75231 Paris Cedex 05, France 

Enrico Cavalli 

Dipartimento di Chimica Generate ed Inorganica, Chimica Analitica e Chimica, Fisica, Universita di Parma, Viale G. P. Usberti 17/a, 43100 

Parma, Italy. 

Marco Bettinelli 

Dipartimento Scientifico e Tecnologico, Univ. Verona Strada Le Grazie 15, 37134 Verona, Italy 



Abstract 

We demonstrate that Zeeman ground-state spin levels in Nd''^:YV04 provides the possibility to create an efficient 
A-system for optical pumping experiments. The branching ratio R in the A-system is measured experimentally via 
absorption spectroscopy and is compared to a theoretical model. We show that R can be tuned by changing the orien- 
tation of the magnetic field. These results are applied to optical pumping experiments, where significant improvement 
is obtained compared to previous experiments in this system. The tunability of the branching ratio in combination 
with its good coherence properties and the high oscillator strength makes Nd^^ :YV04 an interesting candidate for 
various quantum information protocols. 

Keywords: Rare-earth ions, Zeeman transitions, branching ratio, optical pumping 
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1. Introduction 

Rare-earth (RE) ion doped solids have recently 
gained interest in the context of quantum information 
processing applications, particularly as storage devices 
for light at the single photon level jlj]. These quantum 
memories JB |3l |4|, |5l |6t] have a role in future quantum 
networks \ T\, particularly in quantum repeaters |^ 
that can provide a scalable solution for long-distance 
quantum communication. In that context RE ion doped 
solids are interesting due to their long coherence times 
at cryogenic temperatures [10], which could provide co- 
herent storage of single photons on long time scales. 



In quantum memory schemes based on electromag 
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netically induced transparency (EIT) 111 IL 11211 . con 



trolled reversible inhomogenenous broadening (CRIB) 
or atomic frequency combs (AFC) 1 16], the 
system is initialized by preparing the atoms in a well- 
defined ground state (state initialization). This is often 
realized by the usage of a A-system, which consists of 
two ground state levels (typically two spin levels) that 
are optically connected to one excited state. The state 
initialization (or spin polarization) is done by optically 
pumping atoms to the excited state via one leg of the 
A-system and then waiting for the atoms to de-excite to 
the ground state via the other leg of the A-system. The 
capability to perform good spin polarization, with close 
to 100% of the spins being transferred to one of the spin 
states, is crucial to all quantum memory schemes. The 
efficiency of the spin polarization depends on several 
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parameters; (1) on the ground state and excited state 
lifetimes and (2) on the branching ratio R, which de- 
scribes the probabiUty of an excited ion to change spin 
level while decaying. Good branching ratios and long- 
lived ground states can be found in non-Kramers ions as 
3raseodymium [17,, 18.1 . europium | IQ^^QJ^ and thulium 



prase 

m 



2311 ■ It has been considerably more difficult to 



obtain similar results in Kramers ions such as erbium 



|24|,|25t] and neodymium |26], due to short-lived ground 
states and poor branching ratios. Yet, these Kramers 
ions are still interesting as candidates for quantum mem- 
ories due to their long optical coherence times ||27|1 and 
practical wavelengths where diode lasers can be used. 

Here we demonstrate that Zeeman ground-state spin 
levels in Nd^^ :YV04 offers a magnetically tunable 
branching ratio R, resulting in an efficient A-system. 
The branching ratio is found experimentally by measur- 
ing the inhomogeneous absorption profile in a moder- 
ate magnetic field of 0.3 T, where the relative probabili- 
ties of the Zeeman transitions can be directly measured 
from the inhomogeneously broadened absorption spec- 
trum. The branching ratio is measured for two magnetic 
field orientations and the results are compared to calcu- 
lations based on crystal field theory and an effective spin 
§ =1/2 model. Using the found branching value we 
experimentally demonstrate efficient spin initialization 
(97.5% spin polarization) by optical pumping in the A- 
system. These results are interesting in view of the good 
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optical coherence times found in Nd^^ :YV04 12 
and the recent demonstration of light storage at the sin- 
gle photon level in this system iQl- 

2. Theory 

We here investigate the '^19/2 —> '^¥3/2 transition at 879 
nm of Nd^^ ions doped into a YVO4 crystal. This crys- 
tal is uniaxial, and Nd^^ ions substitute for Y^^ ions in 
sites of D2d point symmetry. In a crystal field of this 
symmetry, the ground state '^19/2 splits into five Kramers 
doublets (Zi to Z5) and the ^F3/2 excited state splits into 
two Kramers doublets (Yi and Y2). The transition we 
consider here is between Zi and Yi, see Fig. [T] Un- 
der an applied magnetic field the ground and excited 
state doublets split into two spin sub-levels via a first- 
order Zeeman interaction. This gives rise to a four-level 
system with two ground levels and two excited levels, 
and the associated four possible transitions hereafter de- 
noted a,b,c and d (see Fig. [T)- 

The calculation of the branching ratios between the 
different transitions is based on the electronic wavefunc- 
tions of Nd^^ ions in YVO4 obtained from the diago- 
naUzation of the free-ion and crystal-field Hamiltoni- 



Y, 



abed 



Bo = 



i) 

4'3> = l*-3,2> 
i>,> = \'i'.,n> 



B„//c 



ii) 

l|)4> =-5*|T+3,,>+ Y*|^.3,2> 
*3> = y|^.3/2>+ 5|'P.3,2> 

(^2 > = - P* I T,„,> + a* I T „2> 
*1> = «l^+,/2> + P|^-l/2> 

Bo 5^0 



Figure 1: Investigated transitions, levels and wavefunctions of Nd'^ 
for a magnetic field Bq along the c-axis (i) or at an arbitrary angle 8 
(ii). 



ans, using the parameters obtained by Guillot-Noel et 
al. in Ref. 12811 . This calculation shows that for the ''I9/2 
(Zi) and for the *F3/2 (Yi) Kramers' doublets, the elec- 
tronic wavefunctions are of the form (neglecting term 
and 7-mixing): 



%/2, M = +i) + ^2 



+^3 



79/2, M = +- 



'I9/2, M ■■ 



(1) 



and 



'3/2 



M : 



(2) 



where M is the projection of the total J moment and the 
Aj are complex coefficients. In a D2rf site symmetry, for 
L J multiplet, the wavefunctions are a sum of states 



, 2S + I 
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^^Lj,fi + 4-rrij where fi is the crystal field quantum 
number |29] and m an integer such that -J < fi+4m < J. 
For the '*l9/2 (Zi), we have = +1/2 and for the '*F3/2 
(Yi), fi - ±3/2 as can be seen from Eqs. [T]and|2] With- 
out an external magnetic field, the states characterized 
by +fi and -fi are degenerate because of time reversal 
symmetry (Kramers' theorem). 

Under an external magnetic field Bq, the electronic 
Zeeman interaction Hamiltonian Hez has to be added to 
the free ion an d cry stal field Hamiltonians, where Hez 
is expressed as i3(^\ 

n 

Hez = /JBo- ^ (1/ + gsh) = jSBq- (L + g,S) . (3) 
Here fi is the electronic Bohr magneton, 1, and s, are the 
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individual orbital and spin momentum of the electrons 

n n 

and L = 2 Ij et S = 2 S;- gs is the gyromagnetic ratio 

of the electron spin which is equal to 2.0023. Under an 
external magnetic field, each Kramers' doublet split into 
two electronic singlets. 

From the crystal field calculations, when Bo is paral- 
lel to the c-axis, the wavefunctions of the splitted levels 

are given by and j*!*-^) for the lowest and the 

highest level of the (Zj) doublet and by ^^-I'j and 

for the lowest and the highest level of the '*F3/2 
(Yi) doublet (Fig. [T]i). A general orientation of the ex- 
ternal magnetic field Bq, mixes the and l^-fi"^ states 
and we have for the "^19/2 (Zi) doublet two splitted levels 
\(pi) (low energy) and \(p2} (high energy) given by (Fig. 
fflii): 



Ai) = a|*i'+.)+yS|>F_i) 



(4) 



and for the '^F3/2 (Yi) doublet, two splitted levels \<pi) 
(low energy) and 1^4) (high energy) given by: 



103) = r 1^+3) + 



(5) 



In the following, two branching ratios wiU be discussed: 
^^11 = (6) 



and 



Ri. 




(7) 
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where P\\ is the transition dipole operator along the c- 
axis and P± is the operator perpendicular to the c-axis. 
If we neglect /-mixing, the AJ = 0, +1 selection rule 
shows that the "^19/2 (Zi)^"*F3/2 (Yi) transition has no 
magnetic dipole component. Thus, f y or P± is an elec- 
tric dipole operator corresponding to an electric field po- 
larized parallel (Ey) or perpendicular (Ej^) to the c-axis. 
The weak /-mixing appearing in crystal field calcula- 
tions could explain some experimental results (see Sec- 
tion |4T1). However, this effect is small and will not be 
taken into account in the following. The selection rules 
which hold for electric dipole transitions between the 



states are: 



unless fi - 1^' is evei^S) 
unless jJ^- is odd.(9) 



Moreover, the matrix elements for Py or f ^ are related 
by: 

{^,\ |^_^-) = - {^^,\ 1^/.') (10) 

(^;.|?tuK') = (^-,|/^u|^-,'). (11) 

These relations can be deduced from the symmetry 
properties of the electric dipole operator with respect 
to the crystal field and the time reversal operator 13111 . 

The branching ratios defined in Eqs. |6]and|7]are then 
equal to: 



and 



\a*6 + /3*yf- 
\a*y-(5*6\' 



(12) 



(13) 



The a, /3, y, 6 coefficients are determined by the ori- 
entation of the external magnetic field Bq. They can be 
calculated directly by diagonalization of the sum of the 
free-ion, crystal-field and electronic Zeeman Hamilto- 
nians or by using an eff'ective spin-Hamiltonian model. 
The latter gives more accurate results when the exper- 
imental ^-factor values are known, as is the case here. 
Mehta et al. |32] measured a linear Zeeman eff'ect for 
magnetic fields up to 6 T. This implies that at moder- 
ate magnetic field, as the one we use in our experiments 
(0.3 Tesla), each Kramers' doublet Zi and Yi can be de- 
scribed by an effective spin S- 1/2. Quadratic splitting 
due to the mixing of excited Kramers' doublets with Zi 
and Yi are three order of magnitude smaller than the lin- 
ear terms. We can thus perform a first order calculation 
to determine a, /3, y, 6 coefficients by diagonalization of 
the following spin-Hamiltonian: 122113211 



(14) 



where g is the ^-factor of the doublet under study and 
S the spin 1/2 operator Each ground and excited dou- 
blet Zi and Yi is characterized by the principal values 
g_i_ and g\\ of the ^-factor since the latter is axial for a 
D21/ symmetry. We used the following experimental val- 
ues 4l9/2(Zi):|^^| = 2.361, = 0.915 Q; ^Fs/iiYi): 
\gj_\ = 0.28 (2^, \g\\\ = 1.13 P?]. To be consistent with 
Eqs. |4]and|5] g-factor values must be negative. This 
implies that the = -1/2 spin level is equivalent to 
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Figure 2: Theoretical Sy branching ratio as a function of the angle 6 
between the c-axis and the magnetic field. The light is polarized along 
the c-axis (0=0°). 

the ) State and to the |*^'+|) state for the '^19/2 (Zi) 
and '^Fs/a (Yi) doublet, respectively. 

In the experiments described below, Bq lies in the 
plane containing the c-axis and the E electric field. In 
this case, all a, /3, y, 6 coefficients are real, so that 
^11 - A. plot of /?|| as a function of the magnetic 

field orientation is presented in Fig. |2] At 6=0° (i.e. 
Bo along the c-axis), Nd^^ wavefunctions correspond 
to pure 1^^^ states (Fig. [Tli). In this case, the selec- 
tion rules (Eqs. [8]l show that the a and d transitions are 
forbidden for light polarized along the c-axis. Accord- 
ingly, the R\\ ratio vanishes. On the other hand, when 
Bo is directed away from the c-axis, Nd^^ wavefunc- 
tions are mixture of \^'f,'^ states and all transitions are 
allowed (Fig. [T]ii) until 6 = 90° where the wavefunc- 
tions are again pure j*!*^,^ states. is predicted to reach 
its maximum value, /?|| =0.278, for 0=51° (seeFig.|2]i. In 
the experiments presented in the following sections the 
configuration 0=45° was used, where T^y =0.270, close 
to the maximum value. 

3. Experiment 

The YVO4 single crystals doped with Nd^^ were 
grow n by spontaneous nucleation from a Pb2V207 flux 
il. Reagent grade PbO and V2O5, Nd203 (99.99 %) 
and Y2O3 (99.99 %) were used as starting materials in 
suitable amounts. The doping concentration was 0.001 
% (Nd/Y nominal molar ratio). The batch was put in a 
50 cm^ covered platinum crucible and heated to 1 300 °C 
inside a horizontal furnace. After a soaking time of 12 
h, the temperature was lowered to 850 °C at a rate of 3-4 



°C/h, then the crucible was drawn out from the furnace 
and quickly inverted to separate the flux from the crys- 
tals grown at the bottom of the crucible. The flux was 
dissolved by using hot diluted nitric acid. YVO4 crys- 
tallizes in the 141/amd space group, with cell parameters 
a=b=7.118 Aandc=6.289AandZ=4 |.34ll. 

The crystal sample was cooled to a temperature of 
about 2.8 K in a pulse-tube cooler A permanent mag- 
netic field of 310 mT was applied along different di- 
rections. Light at a wavelength of 879 nm excited the 
"^19/2 —> '*F3/2 transition in the neodymium ions. The 
light pulses were created by gating an cw external cav- 
ity diode laser (Toptica DL 100) with an acousto-optic 
modulator (AOM), which was in a double pass config- 
uration. Pulses of duration 10 ms were created during 
which the laser frequency was scanned over about 15 
GHz by changing the cavity length using a piezo. After 
the AOM the light was split into three beams. A few 
percent were coupled into a Fabry-Perot-cavity for rela- 
tive frequency calibration during the scan. Another few 
percent were used as a reference beam for intensity nor- 
malization, which was detected by a photodiode (Thor- 
labs PDB150A) before the crystal. The third beam was 
focused onto the crystal and was detected in transmis- 
sion by a second, identical photodiode. This signal was 
normalized by the reference signal in order to account 
for intensity changes during the scan. In order to con- 
trol the polarization of the light incident on the crystal, 
a /1/2-plate was installed directly before the cooler 

4. Results and Discussion 

4.1. Measurement of the branching ratio 

The branching ratio R can be determined by measur- 
ing the absorption coefficients of the transitions a,b,c 
and d defined in Fig. [1] We could directly measure these 
by applying a relatively strong magnetic field (0.31 T), 
which splits the Zeeman levels sufficiently so that the 
different transitions could be resolved, despite the 2 
GHz inhomogeneous broadening. We thus measure the 
transmission spectrum while scanning the laser and cal- 
culate the corresponding absorption spectrum. Note that 
we define the absorption depth d in terms of intensity 
attenuation exp(-t/) through the sample. In Fig. |3] we 
show absorption spectra taken with the magnetic field 
parallel 0=0° and at 0=45° to the crystal symmetry c- 
axis. In these measurements the light was polarized par- 
allel to the c-axis for maximum absorption ^2^. 

For the 0=0° orientation, see Fig. |3^, the b and c 
transitions are not resolved, but we can resolve two 
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Frequency (GHz) Frequency (GHz) 

Figure 3: Absorption depth as a function of frequency for a magnetic 
field (a) parallel 8=0° and (b) with an angle of 0=45° to the crystal 
c-axis. In both measurements the light polarization was chosen to 
be parallel to the c-axis. The transitions a,b,c,d are marked with the 
arrows. 



very weak a and d transitions at the expected posi- 
tions. For this orientation we measured a branching ra- 
tio R\\{exp) ~0.05. Note that this value is very approxi- 
mate because of the difficulty of measuring the strengths 
of the weak a and d transitions in Fig. [3^ (particularly 
considering the varying background). The theoretical 
model predicts a branching ratio /?||(,/,)=0 for a magnetic 
field oriented along the crystal c-axis (see Sec. |2]i. The 
small contributions from the a and d transitions that are 
seen in Fig|3^ could be due to a slight misalignment of 
the magnetic field, or due to a small magnetic dipole 
character for the *l()/2(Zi) ^F3/2(Yi) transition which 
does not vanish for this orientation. Indeed, the crys- 
tal field calculation shows that some /-mixing occurs in 
the '*F3/2 excited state with a small contribution of the 
multiplet. The \)/2(Zi) transition with 

A7 = - 1 has a magnetic dipole character. 

For a magnetic field oriented at 0=45° to the crys- 
tal c-axis we can clearly resolve all four transitions, see 
Fig. |3j). Note that there is a slight intensity asym- 
metry in the a/d and b/c transitions, which is not ex- 
pected. This is due to variations in the laser intensity 
during the frequency scan, although most of the varia- 
tion was compensated for using the reference beam (see 
Sec. |3]l. From this spectrum we measured a branching 
ratio of R\\(exp) -0-29 +0.02 (see below on data evalua- 
tion), which is in close agreement with the theoretical 
value 7?||(,/,)=0.27 (see Sec. |2]i for a magnetic field at 
61=45°. 

The theoretical model also predicts that the branching 
ratio is inverted for a polarization orientated perpendic- 
ular to the crystal c-axis, i.e. = i/R±- We inves- 
tigated this by recording spectra while varying the light 
polarization angle. To extract the absorption coefficients 
we fit each spectrum to a simple model consisting of a 
sum of four Voigt functions, with (as free parameters) 
a common absorption depth for the b and c transitions 



dbc, and for the a and d transitions dad- The absorption 
depths as a function of polarization angle are plotted in 

Fig. a 

The dependence of the absorption on the polarization 
angle can be understood by a simple model describing 
the anisotropic absorption. We write the incoming field 
vector (at z=0) as 

EoW = e||£||(0) + e^£^(0) (15) 
= [ey cos(0) + Cj. sin(0)]£'o, 

where ey and are unit vectors, cp is the polarization 
angle with respect to the c-axis, and £0 is the field am- 
plitude strength at z=0. The propagation of the field 
through the crystal is characterized by an anisotropic ab- 
sorption and birefringence. It is important to note that 
the two orthogonal modes transform independently in 
the sample, due to the fact that the linear susceptibility 
is diagonal in the crystallographic axes for a crystal of 
tetragonal symmetry OS] . We can thus write the output 
field amplitude (at z = L) as 

Ei(0) = [ey cos(cf>)e-'«'^ + sm(cf>)e-''-^^]Eo, (16) 

where we have neglected the relative phase factor in- 
duced by the birefringence. This we can do because we 
have no polarization dependent optics after the sample, 
where only the intensity is detected 

Il(.</>) = [cos2(0)e-''« + sin2(<^)e-''-]/o (17) 
:= e-'^^^Io. 

Here we have defined an effective polarization- 
dependent absorption depth d{(p). We fit this formula 
to the experimental data in Fig. |4] (shown as solid 
lines), with d^ and d± as free parameters. As seen 
the agreement is excellent for all polarization angles. 
Note that the weaker a and d transitions show a si- 
nusoidal variation, while the stronger b and c tran- 
sitions deviate from a sinus due to the exponential 
factors in the above formula. The fitted values are 
(i|l„rf=0.75+0.02, d\\bc^2.62+0. 14, d^ad^O.OlO+O.Ol and 
<^±ic=0. 025 +0.01 (all errors are those from the fit shown 
in Fig. |4|i. We can now calculate the branching ratios 
R\\^dnad/diibc^0.29±0.02 and R^^d^ad/d^bc^2.S0±l.2. 
Our measurement confirms the theoretical prediction 
since l/7?j^ =0.36+0. 15 which is equal to R\\ within the 
experimental errors. It should be noted, however, that 
the error in 7?^ is probably larger than above due to 
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Light polarization angle {^) 

Figure 4: Absorption depth d{(p) as a function of light polarization an- 
gle <f) with respect to the crystal c-axis. Squares represent transitions 
b and c {di,c), while circles represent transitions a and d (rfarf). We 
obsei-ve a clear inversion of the transition strength close to the ^=90° 
orientation (see inset), as predicted by the theory (see text). The solid 
lines represent a fitted theoretical model explaining the variation in 
absorption as a function of angle. Note that the magnetic field orien- 
tation is e=45°. 

systematic errors when measuring absorptions at such 
a low level (some percent). Nevertheless, the expect in- 
version in the branching ratio is clearly observed. 

To conclude this section, we have shown that it it pos- 
sible to change the branching ratio by rotating the mag- 
netic field with respect to the c-axis. While such studies 
have been carried out for non-Kramers ions |23], we are 
not aware of similar studies for Kramers ions. Our ex- 
perimentally measured branching ratios at angles 0=0° 
and 45° agree well with the theoretical model presented 
in Sec. |2] A more detailed and quantitative compari- 
son would, however, require more measurements as a 
function of angle. Nevertheless, our work show that 
a theoretical model is useful when choosing promis- 
ing magnetic-field angles. Moreover, we observe ex- 
perimentally the inversion of the branching ratio when 
changing the Ught polarization with respect to the c- 
axis. 

4.2. Optical pumping 

Here we will show that the R value at the magnetic 
field orientation 0=45° makes it possible to perform 
efficient optical pumping. In optical pumping experi- 
ments the branching ratio must have a non-zero value 
(R > 0) in order for the ions to have a certain probabil- 
ity of changing spin state while decaying from the ex- 




Frequency (GHz) 

Figure 5: Absorption spectra before (black line) and after (red line) 
optical pumping. The optical pumping creates a 20 MHz wide trans- 
mission window in the inhomogeneous absorption profile (indicated 
by the arrow). Note that the pumping also creates another hole in 
the leftmost absorption peak, since this transition stalls from the same 
spin level. The pumped ions are moved into the two right absorption 
peaks, where increased absorption is observed after pumping. Inset: 
Zoom of the created transmission window. 

cited state. How large it must be for the optical pumping 
to be efficient depends on the excited state lifetime and 
the ground state population lifetime llzsll . We here show 
that the branching ratio at 0=45° allows for much more 
efficient optical pumping than previous experiments at 
0=0° S. 

The goal was to create a wide and deep transmission 
window within the inhomogeneous absorption profile 
by optical pumping. The total pump duration was 100 
ms during which we made 1000 frequency sweeps over 
20 MHz. The resulting absorption profile was measured 
in the same way as before. The probe pulses were much 
weaker in intensity in order not to modify the absorption 
during the frequency scan. From the absorption profile 
the percentage of ions left in the initial state could be es- 
timated. The orientation of the polarization was 0=20°. 

In Fig. |5]we show the absorption profile measured 
1.3 ms after the burning pulse. This delay allowed for 
all ions to decay back to the ground states (the excited 
state lifetime being ~lOQfj.s). When optical pumping is 
applied (red fine) a 20 MHz frequency window is cre- 
ated in the absorption (no pumping is represented by 
the black line). In Fig. |6]we show the percentage of 
ions left in the ground state, calculated from the absorp- 
tion spectra, as a function of the delay between the end 
of the pump sequence and the time the scan reached the 
position of the bottom of the created transmission win- 
dow. When extrapolating the data to zero delay we find 
that only 5% of the ions are left in the initial state. In 
terms of the total population, which is initially equally 
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Figure 6: Fraction of remaining ions in the pumped spin state as a 
function of the delay time after the optical pumping. The solid line 
is a fit with a time constant of 18 ms. Extrapolating the fit to zero 
delay yields about 4 to 5 %. Inset: Fraction of ions as a function of 
the pump pulse duration. The highest pump efficiency is obtained for 
pulse durations longer than 100 ms. 



distributed over the two spin states, only 2.5% of the 
total population is left in the initial state. This corre- 
sponds to a spin polarization of 97.5%. Note that these 
numbers only apply to the ions that absorb within the 
spectral window created. We can compare this optical 
pumping experiment to the one we made in a previous 
work [26], where we had the magnetic field oriented at 
0°. In that experiment as much as 15-20% of the ions 
were left in the initial state (for a comparable magnetic 
field strength). We can thus conclude that the optimiza- 
tion of the branching factor R unproved the efficiency of 
the optical pumping by a factor 4-5. 

From the time-resolved absorption measurement. 
Fig. |6] we found that about 50% of the initial absorption 
was recovered with a time constant of 1 8 ms. This decay 
time is given by spin population lifetime, since this de- 
cay is much longer than the excited state lifetime of 100 
jt/s. The remaining 50% of the initial absorption recov- 
ered with a significantly slower time constant. Measure- 
ments at delay times of up to 800 ms gave a decay time 
of about 320 ms and at these delays the absorption tends 
towards the original absorption. It seems likely that this 
process is due to population trapping in hyperfine states 
of the isotopes '"^^Nd and ''♦^Nd (12% and 8% natural 
abundances, respectively). An isotopically pure crystal 
would, however, be required in order to study this mech- 
anism for population trapping. Further improvement of 
the optical pumping process should also be possible by 
higher cycling rate through stimulated emission on an 



auxiliary transition ll25ll . Stimulated emission on for ex- 
ample the well-known 1064 nm laser transition in Nd 
should decrease the excited state lifetime considerably. 
In addition one could force spin-mixing of the excited 
state Zeeman levels via radio-frequency excitation to ar- 
tificially increase the branching ratio 12511 . 

5. Conclusions 

We have shown that Nd^^ :YV04 can provide an ef- 
ficient Zeeman-level A-system by tuning the magnetic 
field orientation. The branching ratio is extracted by 
measuring the absorption in a moderate magnetic field 
where the inhomogeneously broadened Zeeman tran- 
sitions can be spectrally resolved. We found that the 
branching ratio can be tuned by optimization of the an- 
gle between the external magnetic field and the crystal 
symmetry c-axis. We find good agreement between a 
theoretical model and our experimental data. 
Finally we used our findings about the branching ratio in 
order to improve spin preparation experiments via opti- 
cal pumping in this material. We created a 20 MHz wide 
transmission window, where only 5 % of the ions were 
left in the initial spin state. Further improvement of the 
optical pumping rate should be possible with spin mix- 
ing and stimulated emission techniques. 
In addition to the findings presented here it is known 
that neodymium shows one of the highest absorption 
of all rare-earth ions, a short excited state lifetime of 
about 100 fis [i22i1 and good optical coherence properties 
iEtI E6ll . Furthermore stable diode lasers and efficient 
single photon counters are available at the transition of 
interest in neodymium, which makes this system inter- 
esting for various quantum information and processing 
protocols. 
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